function  y = buena(a,b,c)

  z = sqrt(b^2 -4*a*c);
  format long 
  x1 = (-b + z)/(2*a) * (b - z)/(b -z)
  x2 = (-b - z)/(2*a) * (b - z)/(b -z)

endfunction


